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Abstract 

It is known that exp[iX{QiPi — i/2)] is a unitary single-mode squeezing operator, where 
Qi,Pi axe the coordinate and momentum operators, respectively. In this paper we employ 

n 

Dirac's coordinate representation to prove that the exponential operator 5„ = exp[iA ^ {QiPi+i + 

Qi+iPi))], (Qn+1 = Qi, Pn+1 ~ Pi), is a n-mode squeezing operator which enhances the stan- 
dard squeezing. By virtue of the technique of integration within an ordered product of operators 
we derive Sn's normally ordered expansion and obtain new n-mode squeezed vacuum states, its 
Wigner function is calculated by using the Weyl ordering invariance under similar transforma- 
tions. 
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1 Introduction 

Squeezed state has been a hot topic in quantum optics since Stoler [1] put forward the concept of the 
optical squeezing in 1970's. 5*1 — exp [iA {QiPi — i/2)] is a unitary single-mode squeezing operator, 
where Qi, Pi are the coordinate and momentum operators, respectively, A is a squeezing parameter. 

Their variances in the squeezed state |0) =sech^/^Aexp — ^aJ'^tanhA |0) are 

AQi = \e'\ APi = \e-^\ (AQi)(APi) = \. 

Some generalized squeezed state have been proposed since then. Among them the two-mode squeezed 
state not only exhibits squeezing, but also quantum entanglement between the idle-mode and the 
signal-mode in frequency domain, therefore is a typical entangled states of continuous variable. In 
recent years, various entangled states have attracted considerable attention and interests of physists 
because of their potential uses in quantum communication |2] . Theoretically, the two- mode squeezed 
state is constructed by acting the two-mode squeezing operator S2 — exp[A(aia2 — aja^)] on the 
two- mode vacuum state 1 00) [31 SI [S], 

S2 1 00) = sechA exp \-a\al tanh aI |00) . (1) 
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We also have S2 = exp[iA((5iP2 + Q2P1)] , where Qi and Pi are the coordinate and momentum 
operators related to Bose operators (ai,a|) by 



Q, = {a, + al)/V2, P, = (a, - al)/{V2i) 
In the state S2 |00), the variances of the two- mode quadrature operators of light field, 



(2) 



X=(gi+Q2)/2, ??= (Pi+P2)/2, [x,<p] = -. 



take the standard form, i.e. 



(3) 



1 



1 



(00| SlX''S2 |00) = -e-^-', (00| S']*p-'S'2 |00) = -e^^, and (AX)(A<p) = -. 



(4) 



On the other hand, the two-mode squeezing operator has a neat and natural representation in the 
entangled state \r]) representation [5], 



52 = 



TT/i 



where 



makes up a complete set 



Iry) = exp(-i + Tja\ - r,*4 + a\4) |00) 



(5) 
(6) 



|r/) was constructed according to the idea of quantum entanglement innitiated by Einstein, Podolsky 
and Rosen in their argument that quantum mechanics is incomplete [7]. 

An interesting question naturally arises: is the n-mode exponential operator 



Sn = exp 



iA V(Q,P,+i +0,+iPO 



(On+i = Qi, P„+i=Pi), n^2, 



(7) 



a squeezing operator? If yes, what kind of squeezing for n-mode quadratures of field it can engenders? 
To answer these questions we must know what is the normally ordered expansion of Sn and what is 
the state S'n |0) (|0) is the n-mode vacuum state)? In this work we shall analyse Sn in detail. But how 
to disentangle the exponential of S'n? Since the terms in the set QiPi+i and Qi^iPi (i = 1, 2, • • • , n) 
do not make up a closed Lie algebra, the problem of what is SnS normally ordered form seems 
difficult. Thus we appeal to Dirac's coordinate representation and the technique of integration 
within an ordered product (IWOP) of operators [HI to solve this problem. Our work is arranged 
as follows: firstly we use the IWOP technique to derive the normally ordered expansion of 5„ and 
obtain the explicit form of 5„ |0) ; then we examine the variances of the 71- mode quadrature operators 
in the state Sn |0), we find that Sn causes squeezing which is stronger than the standard squeezing. 
Thus Sn is an n-mode squeezing-enhanced operator. The Wigner function of Sn |0) is calculated by 
using the Weyl ordering invariance under similar transformations. Some examples are discussed in 
the last section. 



2 Normal Product Form of Sn derived by Dirac's coordinate 
representation 



In order to disentangle operator Sn, let A be 



/ 1 
1 1 

1 



1 \ 






(8) 



V 1 ••• 1 / 
then Sn in ([7|) is compactly expressed as 

Sn = exp[±\QiAijPj], 



(9) 



here and henceforth the repeated indices represent Einstein's summation notation. Using the Baker- 
HausdorfF formula, 

e^Be-^ = B + [A, B] + 1 [A, [A, B]] + 1 [A, [A, [A, B]]] + ■■■ , 



3! 



we have 



S-'QkSn ^ Qk~ XQ^A,k + ^iA2 [Q^A,PJ,QlAlk] 



S-^PkSn = Pk+ XAk^P^ + ^IA" [A^Fj, QM/m^r, 



(10) 



(11) 



From Eq. fTOl) we see that when Sn acts on the n-mode coordinate eigenstate \q) , where q = 
{qiiQ2, - ■ ■ ,<ln), it squeezes \q) in this way: 



Sn \q) = |A|i/' |Ag) , A = e-^^, |A| = det A. 
Thus Sn has the representation on the Dirac's coordinate basis (g|pi?] 

Sn = j dC'qSn \q) {q\ = ^^'^ J d^'q |Ag) {q\ , Sl = S'^ 

since J d^q \ q) {(f\ = \. Using the expression oi\q) in Fock space 



(12) 



(13) 



\q) = TT-"/*: 



exp 



10) 



a) = (a{,4, • • • ,4)' 



(14) 



and the normally ordered form of n-mode vacuum projector |0) (0| ~: exp[— fi^a^] : , we can put Sn 
into the normal ordering form, 



Sn = 7r-"/2|A|i/'y"d"g: exp[-ig(l + AA)q + \/2?(Aa^+a) 



- — (aa + a^a^) — a^a] 



(15) 



To perform the integration in Eq. (|15p by virtue of the IWOP technique, using the mathematical 
formula 

'1, 



/<i"„xpMf. + J„]=,"/^,d«f,-./»exp 



-vP-^v 



(16) 



then we derive 



det A 



1/2 



exp 



is* (AiV- 



detiV 

: exp [at (AiV^^ - ^) a] : exp 



A-/) at 



-S(iV-i-/)a 



(17) 



where = (1 + AA)/2. Eq.([T7]) is just the normal product form of Sn- 



3 Squeezing property of Sn |0) 



Operating Sn on the n-mode vacuum state |0) , we obtain the squeezed vacuum state 



„ , , /detA\^/^ 



1, 



St AN-^A~I] at 



(18) 



Now we evaluate the variances of the n-mode quadratures. The quadratures in the n-mode case are 
defined as 



v2n . , 

2 — 1 ^ l — l 

2 / v2\ 



(19) 



obeying [Xi, X2] = |. Their variances are (AX;) = {Xf^^ — (X^) , i = 1,2. Noting the expectation 
values of Xi and X2 in the state Sn |0), (Xi) — (X2) = 0, then using Eqs. (fTO|) and (llip we see that 
the variances are 



n n 



(AXi)^ (0|5-ix2 5„|0) = i-(0|5-i5^Q,X^Q,5„|0 



j=i i=i 



^ (0| ^ Q.(e-^^)fe. 5;](e-^^),,gH0) 



similarly we have 



2J ' 



^,3 



(AX2)' ^ (0| ^-ix|5„ |0) ^ -E [(AA)-i 



(20) 



(21) 



Eqs. (j20p - (j2ip are the quadrature variance formula in the transformed vacuum state acted by the 
operator exY)[±\QiAijPj\. By observing that A in ([9]) is a symmetric matrix, we see 



n 



2^'n, 



1,3 



» 3 



(22) 



then using AA = AA, so AA = e ^(^+-4) ^ ^ symmetric matrix, we have 

E (AA)^ . = E ^^E [(^ + ^ . , = - E ^^2- = ne--. 



(23) 



and 



„4A 



It then follows 



1 " 



,-4A 



e 



4A 



(24) 

(25) 
(26) 



This leads to (AXi)(AX2) — \, which shows that Sn is a correct n-mode squeezing operator for the 
n-mode quadratures in Eg. p^ . Furthermore, Egs. f^ and (^5)1 clearly indicate that the squeezed 
vacuum state Sn |0) may exhibit stronger squeezing (e"'*'^) in one quadrature than that (e~^'^) of 
the usual two-mode squeezed vacuum state. This is a way of enhancing squeezing. 



4 The Wigner function of Sn |0) 



Wigner distribution functions |12j of quantum states are widely studied in quantum statistics and 
quantum optics. Now we derive the expression of the Wigner function of Sn |0) . Here we take a 
new method to do it. Recalling that in Ref. [T3] we have introduced the Weyl ordering form of 
single-mode Wigner operator Ai (gijPi), 



Ai ((7i,pi)= ■(5(gi -Qi)(5(pi -Pi) ■, 



(27) 



its normal ordering form is 



Ai (91, pi) = - : exp 

TT 



(28) 



where the symbols : : and ; ; denote the normal ordering and the Weyl ordering, respectively. Note 
that the order of Bose operators oi and a| within a normally ordered product and a Weyl ordered 
product can be permuted. That is to say, even though [ai,a|] = 1, we can have : aia|: —: a|ai : 
and ; aia| ; — '. a\ai ] . The Weyl ordering has a remarkable property, i.e., the order-invariance of Weyl 
ordered operators under similar transformations, which means 



(29) 



as if the "fence" [did not exist. 

For n-mode case, the Weyl ordering form of the Wigner operator is 

A„ {q,p)^ '-Sfq-Q^sfp-P 



(30) 



where Q = {Qi,Q2,--- ,Qn) and P — {Pi,P2,--- ,-?„)■ Then according to the Weyl ordering 
invariance under similar transformations and Eqs. pop and (llip we have 



Sn^^n {q,P)Sn 



S-'[5[q-Qjs{p-Pj [Sn 

5 (e^^g - q) 6 [e-^^p- p) 
e-q,e-'^p). 



(31) 



Thus using and dSIl) the Wigner function of Sn |0) is 

(0|5-iA„(<f,p)5„|0) 



— (0| : exp[-(e^-^q - Q)' - {er^'^p- P)1 : |0) 
i^exp[-(e^^g)2-(e-^V)'] 



tt" 



■ exp 



(32) 



From Eq.((32|) we see that once the exphcit expression of AA = exp[— A(yl + A)] is deduced, the 
Wigner function of Sn |0) can be calculated. 

5 Some examples of calculating the Wigner function 

For n = 2, form Eq.([7|) we have ^2 = exp [i2A {Q1P2 + Q2P1)] which exhibits clearly the stronger 
squeezing than the usual two- mode squeezing operator S'2. For n = 3, the three-mode operator 

1 1 

S3, from Eq.([n]) we see that the matrix ^ is | 1 1 | , thus we have 

1 1 



U V V 



A A =1 V u V e 

V V u 



2X 



1 1 



,2A 



3e 



4A ' 



3e4^ 3 



(33) 



and I AA) is obtained by replacing A with —A in A A. Thus the squeezing state 5*3 |000) is 



S3 1 000) = A3 exp 



1000), 



where 



Ai = {I - scch2A) tanh A, A2 



sinhSA 



-, A3 = sechA cosh 2 A. 



2 cosh A cosh 2A 

In particular, for the case of the infinite squeezing A — > cxd. Eg. ([55)1 reduces to 



(34) 



(35) 



S3 1 000) ~ exp ■ 



i=l 



1<J 



1 000) 



(36) 



which is just the common eigenvector of the three compatible Jacobian operators in three-body case 
with zero eigenvalues [13], i.e., 



(Pl+P2 + P3)| ),3 =0, iQ3-Q2 



+ ^2 



' S3 

rrii 



0, 



as common eigenvector 

[Pl+P2 + P3,Q3-Q2] =0, 



+ M2Q2 
M3 + Ai2 



nil + '^■2 + ITT-S 



Ql,Pl+P2+P3 



= 0. 



(37) 



(38) 



Since the common eigenvector of three compatible Jacobian operators is an entangled state, the 
state I is also an entangled state. 

By using Eq. p2|) the Wigner function is 

1 2 ^ 
— ^cxp (cosh4A + 2 cosh2A) > Itti 

TT'^ 3 

i— 1 

r I 3 

X exp < - - (sinh 4A - 2 sinh 2A) ^ of 

[ ^ ^=l 

2 ^ 

- - ^ [(cosh4A- cosh2A) aia* + (sinh2A + sinh4A) aittj] + c.c. ^ . (39) 



]>i=l 



For n — A case, the four-mode operator ^4 is 



the matrix A = 



S4 = exp{iA [(Qi + Q3) (^^4 + P2 
, thus we have 



(Q2 + Q4) (Pi + P3)]} 



(40) 



/ 1 1 \ 
10 10 
10 1 

\ 1 1 / 



AA = 



r t s 

t r t s 

s t r t 

\ t s t r J 



(41) 



where r — cosh^ 2A, s = sinh^2A,t = — sinh 2 A cosh 2 A. Then substituting Eq.(f4T|) into Eq. (f32|) we 
obtain 

(01^4-1 A4 {q,p}S4\0) -- 



■ exp < —2 cosh 2A 



^ V (Af + M*) tanh^ 2 A + {R* + R) tanh2A 



.1=1 



(42) 

where M = aia^ + a2a1, R = a\a2 + aia4 + a^aj, + a3a4. This form differs evidently from the 
Wigner function of the direct-product of usual two two-mode squeezed states' Wigner functions. In 
addition, using Eq. (|4T|) we can check Eqs. (p5)) and (p6|) . Further, using Ea. (|4T)) we have 



det iV = cosh^ 2 A. 



Then substituting Eqs 

54 10000) =sech2Aexp 



/ 2 tanh2A tanh2A \ 

tanh2A 2 tanh2A 

tanh2A 2 tanh2A 

\ tanh2A tanh2A 2 / 

into Eq. pT)) yields the four-mode squeezed state [TTl [TJ] , 



\ (4+4) (4 + 4) 



tanh 2A 



10000) , 



(43) 



(44) 



from which one can see that the four-mode squeezed state is not the same as the direct product of 
two two-mode squeezed states in Eq.dl]). 

In summary, by virtue of Dirac's coordinate representation and the IWOP technique: we have 
shown that an n-mode squeezing operator S'„ = exp[iA ^"^j^ (Q^Pi+i-l-Qi+iPi))], (Q„+i = Qii Pn+\ = 
Pi), is an n-mode squeezing operator which enhances the stronger squeezing for the n-mode quadra- 
tures '16]. iSn's normally ordered expansion and new n-mode squeezed vacuum states are obtained. 
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